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Abstract: We describe the implementation of top production and decay processes in the parton- 
level Monte Carlo program MCFM. By treating the top quark as being on-shell, we can factorize the 
amplitudes for top-pair production, s-channel single-top production, and i-channel single-top produc- 
tion into the product of an amplitude for production and an amplitude for decay. In this way we can 
retain all spin correlations. Both the production and the decay amplitudes are calculated consistently 
at next-to- leading order in as- The full dependence on the 6-quark mass is also kept. Phenomcno- 
logical results are presented for various kinematic distributions at the LHC and for the top quark 
forward-backward asymmetry at the Tevatron. 
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1 Introduction 



The advent of the LHC has brought us data samples of top quarks more than an order of magnitude 
bigger than those observed at the Tevatron. Thus it is opportune to re-examine the issue of top quark 
physics to make sure that the theoretical treatment is at a level appropriate to confront this enlarged 
data sample. We have therefore decided to review the implementation of top processes in MCFM [1-4] 
with a view to providing the most sophisticated next-to-leading (NLO) treatment possible within the 
context of the top pole approximation. 

It is now more than twenty five years since the first calculation of the NLO QCD corrections 
to the production of a pair of top quarks at a hadron collider [5]. This calculation, and subsequent 
NLO computations of single top production [6-12], provided predictions for stable top quarks, without 
consideration of their subsequent decay. This issue was addressed in Ref . [2] , where the decay of the 
top quark was included for the case of single top production, incorporating also the effect of NLO 
QCD effects in the decay. For the case of top pair production, predictions including the decay of the 
top quark are now also available [3, 13-15]. Top pair production processes in association with a jet, 
including the decay, have been considered in ref. [16]. In order to consistently include corrections in 
the decay of the top quark it is imperative to also include the effect of NLO on the top quark width, a 
result that has been available for a long time [17]. In addition to these approaches, complete off-shell 
results for the final state produced by leptonic top quark pair decays, i.e. pp — > ve + be~vb, are also 
now available at NLO [18]. 

The purpose of this paper is to report in detail on a formalism to incorporate QCD NLO corrections 
in both production and decay, in a manner which is easily extensible to other top-production processes. 
We shall keep only the diagrams with a resonant top quark propagator, but provide a complete NLO 
treatment of both production and decay stages of the calculation. In particular we will include NLO 
corrections in the decay of the top quark, including from the hadronic decay of the W boson if 
appropriate, and retain for the first time the effects of a non-zero bottom quark mass. The price for 
including the &-mass effects is very modest in most of the calculation, due to the simple structure of 
the matrix elements [19]. Although the top quark is strictly on shell, we allow the consistent inclusion 
of the effects of an off-shell VK-boson in top decay. 

In this paper we will use this formalism to report on updated results for ti pair production, single 
top production with a IF-boson exchanged in the ^-channel, and single top production with a VK-boson 
exchanged in the s-channel. However, as indicated above, the formalism should be easily extensible 
to more complicated processes such as tiH, tiZ and ttW . The lowest order diagrams for top pair 
production and decay are shown in Figure 1 and the lowest order diagrams for single-top production 
and decay are shown in Figure 2. The s-channel process has previously been treated in ref. [2], but 
now includes the 6-quark mass consistently throughout. This paper is the first instance in which NLO 
corrections have been given for the i-channel single top production process in the four-flavour scheme 
with the decay of the top quark included. A four-flavour treatment of i-channcl single top production 
without inclusion of the top quark decay has previously been given in refs. [11, 12]. The top-pair 
production process has the 6-quark mass retained throughout, and in addition includes QCD radiative 
corrections to the IF-boson decay for the case where the top quark decays hadronically. 

The plan of this paper is as follows. In section 2 we review the results from the literature on the 
total width of the top quark, including mass corrections, off-shell W-boson corrections and NLO QCD 
corrections. Section 3 illustrates the method which we use to calculate top production amplitudes, by 
reference to the specific case of s-channel single top production. Section 4 presents the results on the 
amplitudes for top decay, both at the Born level and the NLO level including both real and virtual 
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Figure 2. Lowest order diagrams for single-top production and (semi-leptonic) decay showing f-channel 
production (upper line) and s-channel production (lower line). The double bars indicate that the top quark is 
on shell. 

corrections. In section 5 we present the counterterm that is used to handle the divergences that occur 
in the real gluon contribution to top decay and in section 6 we explain our procedure for ensuring 
that NLO effects are included consistently in top quark decay. Section 7 presents an illustration 
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Table 1. Notation used for the calculation throughout this paper. 



of phcnomcnological consequences for all three top production processes. The appendices contain a 
summary of the spinor notation used throughout this paper and present a complete calculation of the 
top quark width at NLO. 

2 Results for the top quark width 

In the limit of very large top mass (dropping the masses of the W and the b, s, d- quarks), the lowest 
order result for the top quark width is, 

J 

As usual, Gf is the Fermi constant and Vtj are the CKM matrix elements for the t — > j transitions. 
The superscript (oo) indicates that this is the result that would hold in the limit of infinite top mass, 
when the masses of the W^-boson and bottom quark can be neglected. In the following we shall assume 
Vtb = 1, Vi a = Vid = 0. In this approximation the finite mass corrections implicit in the functions / 
and P 3 correct this result as follows, 

To(o J 2 ) = ^ oa) 2P 3 f. (2.2) 

The notation used in this equation and throughout this paper, taken from ref. [20], is laid out in 
Table 1. Many of the quantities in Table 1 are functions of the variable z which is defined as, 

- - {Pt ~ P 2 W? , (2.3) 



■/?). 



where p t and pw are the momenta of the top and the H^-boson respectively. For the specific case 
of the lowest order process we have that z = f3 2 , corresponding to the &-quark being on shell. The 
variables that depend on z, when evaluated for z = /3 2 , are denoted with a bar. Thus we have, 

1 



P 3 ((3 2 )^P 3 = -v/A(l, W 2 ,/3 2 ), (2.4) 

where the function A is also defined in Table 1. We note that in the presence of one-gluon emission 
we have z m > z > /3 2 . 
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The top-quark width is subject to higher order corrections due to gluon radiation, 

r t = r ( w 2 ) + a s r 1 ( w 2 ). 

The result for the correction to the width has been given in ref. [17] and in rcf. [20] , 
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Taking the limit nib — (/? — > 0) we obtain, 

a s r!(w 2 ) = r(°°)^C F {4(l - w 2 ) 2 (l + 2u?) [Li 2 (l - w 2 ) 



2tt 
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-2w 2 (l + w 2 )(l - 2w 2 ) In oj 2 - (1 - w 2 ) 2 (4w 2 + 5) ln(l - lo 2 ) + |(1 - w 2 )(5 + 9w 2 - 6w 4 )| (2.7) 
The effect of the VT^-boson Breit-Wigner propagator on the total width can be included as follows [17], 



du 2 



(l-£w 2 ) 2 +7 2 



r (w 2 ) + a s ri(w 2 ) 



(2.8) 



where £ and 7 are defined in Table 1. The total width for the decay to an on-shell Vy-boson can be 
obtained from Eq. (2.8) by imposing the narrow width approximation, 



r (i) + as r 1 (i) 



As shown in Table 2 the QCD correction results in a relative change in the top quark width, 

« -0.8a s 



(2.9) 



(2-10) 



that lowers the leading order result by about 10%. 

For our purposes we are interested not so much in the total width but rather in the pattern of 
gluon radiation associated with top decay. Therefore we need a more differential rate, including both 
the distributions of the decay products of the W and of the gluon radiation (if present). This will be 
calculated using a point-by-point subtraction technique that is appropriate for such a calculation [21, 
22]. Thus the integrated real and virtual corrections to the decay rate given in ref. [20] are not 



1 We follow closely the presentation in ref. [20]. Eq. (2.6) clarifies the obvious typographical error S — > 8 in Eq. (27) 
of ref. [20]. We thank Andrzej Czarnccki for correspondence on this point. We have slightly modified the form of the 
expression in Eq. (27) of ref. [20] using dilogarithm identities in order to make the finiteness in the /3 — > limit more 
manifest. 
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m w = 80.398 GcV,m b = 4.7 GeV 


m w = 80.398 GcV,to 6 = 


Tfi w [GeV] 
pw iy bw 


1.453518 
-0.7878491 


1.457412 
-0.7972087 


T$ w [GeV] 


1.476596 
-0.7878090 


1.480522 
-0.7971276 



Table 2. Values of the LO top width and the NLO correction to the semi-leptonic width, using Gf = 
1.16639 x 10 -5 and m t = 172.5 GeV. The results are calculated using both the full Breit-Wigner for the width 
(BW) and using the narrow width approximation (NW). In the former case we take Yw = 2.1054 GeV. The 
first column represents the complete result for non-zero W and b masses, while in the subsequent column we 
approximate by setting the b-quark mass to zero. 



sufficient for our purposes. Nevertheless the computation of the radiative corrections to the width 
contains many essential ingredients for our calculation and also provides an important input for our 
result. Therefore in appendix B we will present details of the calculation of the top decay width, in a 
language that will allow generalization to the case in which we are interested. 

There are further known corrections to the top quark width from electroweak effects [23, 24] and 
from NNLO QCD [25]. In this paper we do not require these corrections. Since we treat the top quark 
as strictly on shell, the width of the top quark enters only as a scale factor, to be chosen to ensure 
that we get the correct branching ratio to the appropriate decay channel. We shall comment more on 
this issue in Section 6. 

3 Amplitudes for top production 

Our approach relies on a factorization of the calculation into amplitudes that include the production 
of a top quark and amplitudes that represent its subsequent decay. This method is restricted to the 
case of top quarks that are produced exactly on shell, but includes all spin correlations in the decay 
of the top quark. By factorizing the calculation in this way we have neglected interference effects 
between radiation in the two stages. Since the characteristic time scale for the production of the top 
quark is of order l/m t while the time for the decay is l/r t , in general radiation in the production and 
decay stages is separated by a large time and interference effects are expected to be small, of order 
asTt/mt [26-28]. 

In this section we will illustrate the method using the case of s-channel top production and briefly 
describe the implementation of the production amplitudes for all top production processes that we 
consider. 

3.1 s-channel single top production 

The production of single top by the s-channel process is the simplest to describe. The lowest order 
process is, 

u{-u) + d{-d) -> 6(p c ) + t(p t ) , N 

I (3-1) 

'->• v{v) + e+{e) + b{p b ) 

where the momenta of the particles are shown in brackets in an obvious notation. We consider all 
momenta to be outgoing so that, 

u + d + p c + p t = u + d + p c + v + e + p b =0. (3.2) 
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Figure 3. Lowest order diagram for s-channel single top production. The double bar indicates that the top 
is on mass shell. 



Since our method relies on the assumption that the top quark is exactly on shell, and the bottom 
quarks are also on their mass shells we have that, 

2 2 2 2 2 to o\ 

Pt = m n Pc=Pb = m b- (3-3) 
Thus the amplitude for the process in Figure 3 is, 

M Xb x c =6 idi J ibi / w A d\YM u Xb ( Pb )r lL ^ t+ ^ t) 7p7L V Xc (p c ) / w H7"|e], (3.4) 
2D w {s ud ) im t T t 2D w (Sue) 

and g^r = A\/2m 2 v GF- We have suppressed the colour labels id, i u , % and i c of the d, u, b and b 

quarks on the left-hand side. The matrix is denned by 7l = 1 — = |(1 — 75), the propagator 
factor is, 

D x {s) = (s - ml) + im x T x , (3.5) 

and Sij = (i + j) 2 , Sijk = (i + j + k) 2 where i, j and k stand for four-momenta. The notations |z) and 
\j) denote spinor solutions of the massless Dirac equation. A brief summary of our spinor notation is 
given in Appendix A. Using the identity, 

u\ t {pt)u\ t (p t ) =p ! t + m t , (3.6) 

a, 

we may write Eq. (3.4) as, 

M XbK = J2 ^Tf P ^ < ( 3J ) 

At, it 

where the amplitudes representing the production (P\ t \ c ) and decay (D\ b x t ) of the top quark are 
given by, 

g 2 - g 2 

D\ b \ t = 5i bit u\ b (p b )YlLU\ t {pt) — -j — r(^|7 M |e] = h b i t n Y , u Xb (p b )\v){e\u Xt (p t ) , (3.8) 

g 2 

P\ t \ c = S it iji d i u w , — Ad\Y\u] u Xt (Pth P 7LVx c (Pc) 
ZUwysud) 

= Si t iji d i u 9 J u Xt (pt)\d)[u\vx c (p c ) . (3.9) 
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These equations have been simplified on the right hand side by using the Fierz-like identity, Eq. (A. 5). 
In order to proceed further we shall use the standard trick [29] of decomposing the massive momenta 
into the sum of two massless momenta, p 1 ^ = f + Q-trjf with the constant at given by, 



«t = , | 1 ■ (3.10) 

We may write the massive spinors as combinations of massless spinors as follows, 

u-(pt) = (j^ + mt)I%>7r-T, u +(Pt) = (j4 + m t )\r] t }-J— , (3.11) 

u-(pb) = [r)b\(tfb + m b)-r~u> u+(Pb) = (VbKp'b + m b)-r-rr , (3-12) 
[Vbb\ {iibb) 

v +(Pc) = (i>c - m b )\r/ c )j—^, v-(p c ) = (tfc - m b )\r) c ]j— -j , (3.13) 

V+(Pa) = VlaWa - m t )— !— , V-(p a ) = {rjaWa ~ TUt) ~, : T ■ (3.14) 

[i]aa\ (r}a,a) 

The momenta pt and p a refer to the top and anti-top quark respectively whereas p b and p c refer to 
the bottom and anti-bottom quark. The spin labels of the massless spinors 1 77^ } , 1 77^ ] for i = t,a,b,c 
encode the polarization information of the massive quarks and they are equivalent to helicities only 
in the massless limit. Stripping the overall colour, coupling and propagator factors, 

Px t x c = S iti J idin 9 f V XtK , (3.15) 

J->W\Sud) 

we can evaluate the production amplitudes in Eq. (3.9) to find, 

V __ = - mh ^^ _> 

[VcC] 

V-+ = - (td) [cu] 

(rjt d) [q c u) 

V+- = ~m h m t -. — r -> 

\Vt t) [T]c c\ 

(r) t d)[cu\ (ed) [cu] 

V++ = -m t — — 7 > m t — — — . 3.16 

(r) t t) (t e) 

The final expressions (in which two of the amplitudes vanish) are valid for the special choices, rjt = 
e, f] c = u. Unless otherwise stated, we shall always simplify our amplitudes with these choices of r) t 
and rjc- 

3.1.1 Virtual corrections to s-channel single top production. 

As is apparent from Figure 4 the virtual corrections to the s-channel production are of two types. In 
the lower row we have the vertex corrections to a timelike vector boson (a) coupling to two massless 
fermions and (b) coupling to a massive top quark and a massive bottom quark. These corrections 
are well known and will be considered in more detail in the next section in the context of one-loop 
corrections to the decay of a top quark. Indeed the identical corrections occur (a) in the hadronic 
decay of a W boson, and (b) in the one-loop corrections to the semileptonic decay of a top quark. We 
will therefore defer our discussion of these corrections until later. 
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Figure 4. Diagram for s-channel single top production with virtual gluon radiation. 




Figure 5. Diagram for s-channel single top production with additional gluon radiation. 
3.1.2 Real corrections to s-channel single top production. 

Figure 5 shows the diagrams contributing to s-channel single top production in association with an 
additional radiated gluon. The diagrams in the top row have gluon radiation from the initial state 
quarks, whereas the diagrams in the bottom row have gluon radiation from the final state quarks. The 
two rows of diagrams are separately gauge invariant and do not interfere because of colour conservation. 
In addition to the diagrams in Fig 5 there are also diagrams obtained by crossing the gluon in the top 
row into the initial state. The diagrams obtained by crossing the gluon in the bottom row into the 
initial state are not included because they correspond to the i-channel single-top production process 
and are therefore accounted for elsewhere. We begin by defining the reduced matrix elements for the 
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real radiation from the initial and final states, 

^ilUc = 9w9s (t A ) ldlu D w {s udg )V^ XgXc , (3.17) 
P { x t l g x c = 9w9s (t A ). t . c D w { Sud )vi f X gXo > ( 3 - 18 ) 

where the colour matrices t A are normalized such that t A t B — 8 AB . The spin labels At, A^, A c encode 
the polarization state of the top quark, gluon and bottom anti-quark respectively. After a simple 
calculation we find for radiation from the initial state, 

v (i)__ = o 

v (i) _ {t\d + g\u} [cu] 



V {i)__ = o 



[g u] [g d] 

m t (e\d + g\u] [cu] 
(t e) [g u] [g d] 
(i) _ m b (t d) [g u] 
(gu) [cu] 

l\u + g\c] 
(g u) (g c 

i) m t m b (ed) [gu] 



v (i) (d\u + g\c\ (td) 
'-++ (gu)(gd) 



(te) (gu) [cu] 



v (i) 1JH (d\u + g\c] (e d) 

+++ (te){gu)(gd) 

The results for radiation from the final-state quarks are, 

V U)__ = o 

(/) (td)[cu] f (g\p t \d] (g\ Pc \d} \ (tg)(gd)[cu] 

+ [gd] \(g\ P t\g] (g\ Pc \g]J (g\ P t\g] 
pg_ = 

(/) m t (ed) [cu] ( (g\p t \d] (g\p c \d] \ m t (ge) (gd) [cu] 

+ - + (te)[gd] \(g\ Pt \g] (g[p c \g]) + (te)(g\p t \g] 
V U) _ m b (td) [gu] 2 



[cu] (g\p c \g] 

(/) . ( td ) [cu] ( (u\p t \g] (u\p c \g]\ (td) [eg] [gu] 



(gu) \(g\pt\g] (g\p c \g]J (g\p c \g] 

(/) m t m b (ed) [gu] 2 

(te) [cu] (g\p c \g] 



(/) mt (ed) [cu] f (u\p t \g] (u\p c \g] \ m t (e d) [c g][g u] 

(te)(gu) \(g\ Pt \g] (g\p c \g]J (te)(g\ Pc \g] ' K ' ' 

3.2 Top pair production 

In our current treatment of top pair production we have followed the procedure outlined above and 
expressed the complete amplitude as a product of the amplitude for ti production times the amplitude 
for the decay of the top quark and the anti-quark. In addition we have coded the analytic results for 
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the one- loop amplitudes for top pair production given in the paper of Badger et al., ref. [30]. These 
amplitudes are expressed simply in terms of spinor products. Our previous treatment of the one-loop 
amplitudes [3] was based on the work of Korner et al. [31] in which the results are expressed in terms 
of more complicated spinor strings. We note that the two calculations are in complete numerical 
agreement. We find that using the new expressions of Ref. [30] in MCFM improves the speed of the 
calculation of the virtual corrections by a factor of three. 

3.3 i-channel single top production 

For the ^-channel single top process we use the four-flavour scheme in which the lowest order process 
is, 

q + g ^ q +b + t{-^ ve + b) , (3.21) 

as indicated in the top row of Figure 2. The next-to-leading order corrections to this process, without 
including the top decay products (i.e. qg — > q'bt), were presented in Refs. [11, 12]. In this paper we 
re-use the amplitudes computed for that calculation but extend them to incorporate the top quark 
decay in a fashion similar to the approach presented for the other two cases above. A slight difference 
is that the virtual amplitudes used in Ref. [11, 12] had already been simplified by making a specific 
choice for the vectors that are used to decompose the massive top quark momentum, rj t . In that case 
the standard decay amplitudes that are presented in the following section cannot be used. Rather 
than repeat the calculation of the relevant amplitudes with the canonical choice of this vector used 
throughout the rest of this paper, we have simply used expressions for the Born-level decay amplitudes 
(see Eq. (4.4) in the next section) that are appropriate for this alternative choice of rjt ■ 



4 Amplitudes for top decay 

In this section we list the various decay amplitudes that are used in the calculations presented in 
this paper. Except where otherwise noted, we make the choices rj t — e, rjb — with the choice of rjt 
required in order to match the production amplitudes. 

4.1 Top decay at Born level 

The Born decay amplitudes for the process 

t{ Pt )^b{ Pb ) + V { V )+e+{e) (4.1) 

can be easily evaluated from the expression given in Eq. (3.8). In order to report the results using 
spinor products it is useful to define the massless vectors, t ■ t = b ■ b = 0, using the auxiliary vectors, 
Vb,Vt, 

<* = *•- si*- 

Separating out the colour, coupling and propagator factors, 

D\ b \ t = 5 ibit g ™ — - V Xb x t , (4.3) 
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we find, 

= (bu)[et] 

ml 

(Vbb) 

V ++ = m b m t ^-M^O. (4.4) 
(Vbb) [tilt] 

Making the choices rj t = e (as in the production amplitudes) and rjb = u 1 all but the amplitude D 

vanish. 

As noted above, for the particular case of virtual and real radiation corrections to the i-channel 
production process the amplitudes have already been simplified using r\ t =9, the momentum of the 
gluon in the initial state. One can read off the appropriate decay amplitudes from Eq. (4.4) by making 
the choices rjt — g,rjb = v. We see that in this case the amplitude T> |_ does not vanish. 

The amplitudes for the related charge-conjugate process, 

t(p a ) -> b( Pc ) + e~ (e) + v{9) , (4.5) 

can be obtained from the above amplitudes by symmetry. Denoting the colour-stripped amplitudes 
for this process by T>\ a \ a and introducing the operator C defined by, 

C : t — > a , b — > c , v — > v , e — > e , (ij) •H- [ij] , (4-6) 

then we have, 

V XaK = -C[ZU 6 _ At ] . (4.7) 

Note that the massless vectors a and c that appear in the expressions for the amplitudes are related 
to the massive four-momenta of the top and bottom quarks, p a and p c , by the transformed equivalents 
of Eq. (4.2), 

9 9 

with the choices r\ a = e and r\ c = v. 

4.2 Virtual corrections to top decay 

We now discuss top quark decay with virtual gluon radiation from the t, b line as shown in the top row 
of Figure 6. We give the results for the special choices, r\ t — e, rjb = v for the one loop corrections to 
the decay, 23^ A . The labels At, and At denote the spin state of the bottom and top quark respectively. 
In the massless limit they would correspond to helicitics. The general result for the rcnormalizcd 
one-loop form factor in the case m t ^ 0, rrib ^ is, 

Fti(pb,Pt) = ^u{pbW ^{pb,Pt)u{p t ) , (4.9) 



where r M (pb,p t ) is decomposed as, 

T»{p blVt ) = C F cr(-4V 



^ ^ ^ ^ 



<^0 7 71 + ^0 TlR + ^1 1R + L-l 7L + ^2 1R + C 2 1L 

m t m t m t m t 



+ O(e) (4.10) 
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VZ t VZ b x 





+ 



Figure 6. Diagrams for gluon radiation from the t, fo-line in top decay. Zj, and Zt denote the wave function 
renormalization constants for the external fermion lines. The expression for the Z q ,(q — b,t) is given in the 
appendix, Eq. (B.7). 



The factor cp is the standard prcfactor that appears in dimcnsionally regulated one-loop calculations, 

1 T(l + e)r 2 (l - e) 1 1 



cr 



(4tt) 5 



r(l-2e) 



(4tt) 2 " £ r(l - e] 



(4.11) 



If we ignore the masses of the decay products of the W-bosoia, the C 2 terms will not contribute to 
physical amplitudes. They will be dropped in the following. The results for the remaining coefficients 
are given in Eq. (B.6) of Appendix B.2. The amplitudes , Cf and are finite in the limit e — > 0. 
Removing the common overall factor, 



9w9 2 s ctCf /M 2 V 



V 



we have that , 



(4.12) 



V v __ = <C%{bv) [et} + 



Cf (bv) (eb) [eb] 
~~2 (Ft) 



m t m b C^ 



e v) [v e] m b Cf (b v) \v t] [e b] 



+ 



(et)[vb] m t 2 



[1/6] 



VZ. = { m b C* 



[1/6] 



C[ (vb) [ve] [eb] 



m t C Q 



-L-(vb)(bt)[eb]+m b , 
R (ue) [e b] Cf , . ,_ . , _ C« (i/ e) [e 6] 



vb] [et] 



(vb) [eb] [bt]+m b - 2 



(et) 



^.v I , . , ,, Cf (6z/) [e t>] 2 CP m b , , , 
^ + = Cf <t „> [e6] + ^-^T^ + -f~ <*"> [e6] 



(4.13) 



Because of the vanishing results in Eq. (4.4), T)+\ t = 0, for our particular choice of auxiliary vectors, 



rjt, r\b the one- loop results for T>Y_,T>Y + will not be needed. 
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The virtual amplitudes for the decay of an anti-top quark are related to those for the decay of a 
top quark in the same manner as for the lowest order amplitudes. Namely 

vl aXc = -C[D\_ Xt ] , (4.14) 

where the operation C is defined in Eq. (4.6). 

4.3 Top decay with gluon radiation 

We now discuss top quark decay with gluon radiation from the t, b line as shown in the bottom row 
of Figure 6. In order to report the results using spinor products it is useful to define the massless 
vectors, t ■ t = b ■ b = 0, as in Eq. (4.2), using the auxiliary vectors, e, v, 



t" 



Pb 



2p t .e 2pb.v 
We further define the colour-stripped decay amplitudes through the relation, 



<A 3 A ( =9sf{t A ) ibit 



9w 



■ V R 



(4.15) 



(4.16) 



Dw{s vS ) 

where the colour matrices t A are normalized such that t A t B = 8 AB . The explicit expressions for the 
decay amplitudes are, 



(bu)[te] / (g\ Pb \b] (g\p t \b}\ , (bg)(gu)[te] 



[bg] 

(bv) [te] 



V 



V 



R 



++ 



-+ 



+++ 



(g\pb\g] (g\pt\g) i 

[g\Pb\b) [g\ Pt \by 
(gb) y\g\pb\g) [g\pt\g) ^ 

m b (v g) 2 [et] 

{vb) (g\pb\g\ 
o, 
o, 

m t {vb) [eg] 2 
[et](g\pt\g] ' 
0, 
. 



(g\pb\g] 

M [gt]\ge] 
[g\pt\g) 



(4.17) 



and V^ + _ 



As expected the soft gluon singularities are all contained in the eikonal terms for T>^_ 
These can be simplified using the Schouten identity, Eq.(A.6), to remove the unphysical singularity. 
For example, the amplitude T> contains the singular term, 



1 f(g\Pb\b] (g\ P t\b] 



(g\Mt\g) 



(4.18) 



[bg] \(g\pb\g] (g\pt\g]J (g\pb\g] (g\pt\g] 

Squaring the part of the full decay amplitude that contains the eikonal factor we obtain the standard 
expression, 



A.A\ 



g z s V 2e Tr(t A t 



ig\Mt\g) 



(g\pb\g\ (g\pt\g] 



i 



g 2 ^C F x 



2p t -Pb 



(Pb-g) 2 Pb-gpt-g (pb-g) 2 



\D-r > 

(4.19) 
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and an equal result from the other gluon polarization, D^, . 

This eikonal form is the basis for the counter-term that we use to subtract the singularities resulting 
from radiation of a soft gluon in the top quark decay. We make use of the factorization, 

g 2 s C F ^ 2e ^2 P^g^tiPuPw ,Pb,P g )\ 2 -> X] \ v >-b*t\ 2 (Pt,Pw,Pb) * S(p t .p g ,Pb-Pg,m^ml,pw) , 

(4.20) 

where the lowest order matrix element T> is evaluated for momenta pw, Pb of the M^-boson and 6-quark, 
as described in more detail in Section 5. In practice, we shall not use the complete eikonal factor as a 
counterterm, but rather a simpler expression that differs from the above by finite terms. Using 



(Pt ~Pb- g f = Pw 



and dropping non-singular terms we obtain, 

S{Pt-Pg,Pb-Pg,rnt,ml,Pw) = g 2 ([J, 2 ) e C F x 



Pw 



(pb-g) 2 



Pb-gpt-g 



{pb-gf 



(4.21) 



(4.22) 



The amplitudes for the emission of a gluon from the t-b line in anti-top quark decay are related 
to the amplitudes presented above, in a similar manner to the lowest order case. We have, 



V 



A a A s A c 



c 



V 1 



(4.23) 



where the operation C is defined in Eq. (4.6). We note that, since these amplitudes contain an addi- 
tional gluon, this relation is different by an overall sign from those presented in Eqs. (4.7) and (4.14). 

4.4 Top decay with virtual gluon radiation from the decay products of the W 

If the W boson decays hadronically, 



t(p t ) -> b(p b ) + q{q) + q(q) 



(4.24) 



then we should also include the next-to- leading order corrections to that decay process, as depicted in 
Figure 7. The result for the virtual corrections to the decay of a vector boson into a pair of massless 
fermions is well known [32]. We can write the form factor as, 



where the expansion is very simple compared to the general case (c.f. Eq. (4.10)), 



T"(q,q)=^ L [l + g 2 C FCr C L 
Stripping overall colour-factors as usual, 



0(e) 



jjVw 
Ab A+ 



- _ A g w 9 s c rC F / 

\ t °ibh n f \ \~ 
Dw(Sqq) V - 

all of the amplitudes vanish except one, 



Afa At 



(4.25) 



(4.26) 



(4.27) 



v Vw = C L (bq)[qt] . 



(4.28) 
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+ 




+ 




Figure 7. Radiative corrections due to the interaction of the decay products of the W in top decay. The 
upper line shows the virtual corrections and the lower line the real corrections. 



Note that, for the labelling of momenta indicated in Eq. (4.24), the massless momenta t and b are 
defined by choosing rj t = q and rjb = q in Eq. (4.2). In the four-dimensional helicity scheme the 
coefficient C L takes the familiar form, 



c L = -- 2 --- 

The amplitudes for the related charge-conjugate process, 



7 . 



t(p a ) -> b{p c ) + q{q) + q(q) , 



(4.29) 



(4.30) 



can be obtained by a symmetry operation that differs from that in Eq. (4.6) only by a relabelling of 
momenta. We have, 



-c 



w 



— \t> —At 



(4.31) 



where Cw is defined by, 

C w ■ t->a, b^c, q^q, q^q, (4.32) 
The massless momenta a and c are defined by the choices r) a = q and r] c = q in Eq. (4.8). 

4.5 Top decay with real gluon radiation from the decay products of the W 

The diagrams for real radiation of a gluon in the W decay are shown in the bottom row of Figure 7. 
After removing overall factors, 



jjRw 



9w 



D w (s qqg ) Wt ' 



(4.33) 
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a simple calculation yields the amplitudes, 



V 



v Rw_ 

T<Rw 

u ++- 



V 1 



= m b - 

= 0, 

= o , 



[qt] (b\pw\q] 
[qg] [qg] 

(qb) (q\p w \t] 

(qg) (qg) 
(qg) [<jt] 



q b) [q g] 



u -++ 



-m t - 



( qb) [qg] 
qg) [qt] 



+-+ 

Rw 

++ + 



(4.34) 



where pw = Pq + Pq + Pg- The counterterm in this case is the standard Catani-Seymour dipolc 
associated with final-final radiation from massless partons [22]. 

The equivalent amplitudes for the anti-top quark decay process are given by, 



V 



V 1 



(4.35) 



where the operation Cw is defined in Eq. (4.32). Once again this relation has the opposite sign from 
that in Eq. (4.31). 



5 Counter-term for real radiation 

In this section we describe the implementation of the counter-term that is used to cancel the soft 
singularity resulting from real radiation from the t-b line in the decay of the top quark. The counter- 
term has already been introduced in Eq. (4.20), which makes clear that the lowest order matrix element 
is evaluated at momenta pw and pi, of the W boson and bottom quark respectively. These momenta 
are related to those for which the real radiation matrix element is evaluated (pw, Pb, Pg) as follows. 

The transformed momenta pw and pb are generated by a Lorcntz boost along the direction of 
the W in the top rest frame, following the method of ref. [2] (Section IV). The transformation is 
generalized slightly in order to incorporate the effect of a non-zero bottom quark mass. Thus we still 
have p w = p w with, 



+ 



sinh(a;) 



V(Pt - Pw) 2 -PwPt 
cosh(x) — 1 



PtPw ~ PwPt ) 



(Pt - Pw) 2 -PwPt 



(pt ■ pw (PtPw + PwPt) - Pw PtPt - Pt PwPw) (5- 1 ) 



but pb is now constrained by p 2 — (p t — Pw) 2 = rn\. This constraint determines the coefficients that 
appear in the transformation, 



sinh(x) 
cosh(x) 



1 



2 P 2 tP 2 w 
1 

2p 2 tP 2 w 



- V Hpf,Pw>Pb)pt 'Pw + (Pt +Pw -Pb)y (pt - pw) 2 -p w Pt 



+ (Pt +Pw~Pb)Pt -Pw - JHphPw'Pb) V (Pt -Pw) 2 -PwPt 



(5.2) 
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This Lorentz transformation is also used to determine the modified momenta of the W decay products. 

With the subtraction fully specified, we now turn to the issue of performing the integration of 
this counter-term analytically in order that the poles can be extracted and cancelled against those 
appearing in the virtual calculation. From Eq. (B.14) we may write the phase space for the decay of 
an on-shell top quark as 

d<^ ( ' i \p t - 1 p w ,p b ,p g ) = d$ {2} (p t ;pw,Pb) x [dg(p t ,pw,z)] 

= d^ 2) (p t ;pw,Pb) x [dg(p t ,pw,z)] ■ (5-3) 

The equivalence in Eq. (5.3) follows from Eq. (B.l) in the appendix because d n ~ 1 uj w = d n ^ 1 Lu 1 x, since 
Pw and pw are related by a boost. 

The counter-term may now be integrated in exactly the same manner as the real radiation matrix 
elements that enter the calculation of the total width, as detailed in Appendix B. The integrated 
counter-term corresponds to the sum of the integrals Si , S2 and S3 , expressions for which are given in 
Eqs. (B.23), (B.24) and (B.25). We thus arrive at the final expression for the integrated counter-term, 

[dg(pt,pw,z)] S{p t .p g .,Pb-Pg,m 1 t ,m 1 h ,p\ v ) = Si + S 2 + S3 = 

2 P2 \ 1 



2 \ e 

2 9 2 ctCf 1 



H^)} (i-^) +4+ A( (i - w,)y ' +(i - /3B)p -) 



p 3 



2Y P - 6Y p 2 + 4Y W ln(/3) - 6Li 2 ( 1 - t- ) - 2 Li 2 (1 - P+) + 2 Li 2 (1 - P_) } . (5.4) 



6 Consistent treatment of top quark decay in perturbation theory 

Since we are everywhere treating the top quark as on shell, the full cross section integrated over the 
decay products of the top will be given by the production cross section multiplied by the branching 
fraction to the chosen decay channel. This statement should hold at all orders in perturbation theory 
For example, we can write the differential NLO cross section for single top production followed by 
semi-lcptonic decay schematically as, 

a NLO ( PP -+ ue+b) +X) = K + as a x ) x gfl+^gf , (6 . 1} 

1 + a s i 1 

where 00, Tg\ Tg arc the lowest order contributions to the production rate, semileptonic decay width 
r(i — > i>e + b), and total top width and ageri, ots±i and asFi the corresponding NLO corrections. 
The total top width is given by the sum of the partial widths to the various decay channels, 

T t = J2 T (t^fifjb + X), (6.2) 

i,3 

which factorizes in the narrow width approximation for the W-boson to, 

r t -> r{t^ bw + x) x V r( - w T ^ ^ = r(t -^bw + x). (6.3) 

The quantity T(t — > bW + X) has the perturbative expansion, T(t — > bW + X) = T a + agTi, where 
To, Ti are as given in Eqs. (2.2), (2.6). This is the expression for the total width used in the denominator 
of Eq. (6.1). 
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In order to include the effect of radiation in the decay of the top quark properly, we must be careful 
to ensure that we perform the perturbative expansion in the strong coupling in a consistent manner. 
As written in Eq. (6.1) the NLO calculation includes a contribution of relative order a s corresponding 
to corrections in both production and decay stages simultaneously. We follow the treatment in ref. [14] 

4 



and simply expand Eq.(6.1) in as, discarding terms of order a s or higher, 



(/) 



„NLO „ „ ^0 

a — (Tq x — — - 



+ (TQ X 



(/) 



+ as(T\ x 



dT 



(0 



- asCTo x 



dT ( l) Ti 



(6.4) 



Performing the full integration over the final state, i.e. making the substitutions dT^' — » Iq and 
dTj — > , we recover the NLO cross section for corrections in production only, 



a NLO = ( a(j + as(Ti ) x Br ^ w _^ ue +^ + asao Br(W -> ue + ) 
= (cr + u s <Ji) x Br(W ve + ) 



r (0 
1 1 

1 o 



To 



where the equality follows, for example in the narrow width approximation, because, 

r(0 



ri(t-> bW + X)Br{W 
~ T (t ->• bW + X)Br(W 



(6.5) 



(6.6) 



We note that while Eq. (6.5) is a desirable outcome, it is only true in the case that all degrees of 
freedom associated with the decay are completely integrated out. In the presence of experimental 
cuts, for instance on the leptons or 6-quarks present in the decay, the two calculations (i.e. with and 
without radiation in decay), can predict different cross sections. The above discussion pertains to 
the case of NLO corrections to the production and leptonic decay of a single top quark. For top pair 
production, with both top quarks decaying lcptonically, the last term that appears in Eq. (6.4) appears 
once for each quark so that this correction factor is doubled. 

We note that since the top quark is produced exactly on-shcll, the value of the top quark width is 
important only to ensure that we obtain the correct branching ratio when including the decay. This 
can clearly be achieved by using Tq as shown in Eq. (6.4). This is therefore the value of the width 
used in our code. 

The argument given above also holds, mutatis mutandis, for hadronic instead of leptonic decays 
of the VF-boson. 



7 Phenomenology 
7.1 Input parameters 

The results presented in this paper are obtained with the latest version of the MCFM code (v6.2). 
The electrowcak parameters that we regard as inputs are, 

M w = 80.398 GeV , T w = 2.1054 GeV , (7.1) 
G F = 1.16639 x 10~ 5 GeV" 2 . (7.2) 

The top and bottom quarks have the pole masses, 



m = 172.5 GeV , m b = 4.7 GeV 



(7.3) 



Treatment of V^-boson and &-quark 


olo [fb] 


erjvLO (prod.) [fb] 


(Tjvlo (prod.+dccay) [fb] 


Narrow width, m& = 
Narrow width, m& = 4.7 GeV 
Breit-Wigner, m& = 4.7 GeV 


30.98(2) 
30.78(2) 
30.77(2) 


48.84(3) 
48.61(3) 
48.61(3) 


48.82(3) 
48.60(3) 
48.59(3) 



Table 3. Cross sections in femtobarns for the s-channel single top process, computed with various levels of 
sophistication. No cuts have been applied to the final state. The numerical integration errors are shown in 
brackets. 



For the parton distribution functions (pdfs) we use the sets of Martin, Stirling, Thorne and 
Watt [33]. For the calculation of the LO results presented here we employ the corresponding LO pdf 
fit, with 1-loop running of the strong coupling and a s (Mz) = 0.13939. Similarly, at NLO we use the 
NLO pdf fit, with a s (M z ) = 0.12018 and 2-loop running. 

In our calculations of the top pair and s-channel single top processes we set the factorization and 
rcnormalization scales equal to the top quark mass. For the ^-channel process we employ two scales, 
one for evaluating contributions associated with the light quark line (/i;) and the other for the heavy 
quark line (/i/j). We set [i\ = m t /2 and fj,h = "i t /4, as advocated in Ref. [12]. 

7.2 s-channel single top at the Tevatron 

In order to illustrate some of the features of our calculation, we begin by considering the s-channel 
single top process. Although the thrust of this paper is the inclusion of the top quark decay allowing 
realistic experimental cuts to be applied to the top decay products, it it instructive to first consider 
the case in which no cuts are applied. 

The predicted cross section for the process pp — > t(— > ve + b)b at the Tevatron, computed at 
various levels of sophistication, is shown in Table 3. In this table we present results of performing the 
calculation at LO, at NLO including radiation in production only, and at NLO including radiation in 
both production and decay. Moreover, we consider three different kinematic approximations: setting 
m& = and using the narrow-width approximation for the VF-boson (equivalent to the approach in 
Ref. [2]), using m& = 4.7 GeV but still using the narrow-width approach, and finally using m& = 
4.7 GeV and implementing the full Breit-Wigner form of the W propagator. We first observe that, by 
construction, the two columns of calculated NLO cross sections are in perfect agreement, thanks to the 
presence of the last term in Eq. (6.4). Moreover, the slight adjustment to the width of the VT-boson 
when using the Breit-Wigner rather than the narrow- width approximation, as indicated in Table 2, 
ensures that the final two rows of Table 3 are also in excellent agreement. Finally, we observe that the 
effect of including the 6-mass in this process is very small, resulting in a decreases in the cross-section 
by approximately 0.5%. We note that this effect is not the result of including the mass of the 6-quark 
in the top quark decay but is instead due to the treatment of the b quark. 

Wc now perform the same analysis using a realistic set of experimental cuts. In particular we use 
the cuts employed in a recent search for the Higgs boson using the WH (— > bb associated production 
channel, for which s-channel single top production is an irreducible background [34]. All jets (both 
light and heavy flavour) are clustered according to the anti-/cr algorithm with a distance parameter 
D = 0.7. Wc require that the algorithm finds at least two jets that satisfy, 

p T (jct) > 20 GeV , y(jet) < 2 . (7.4) 

and the charged lepton acceptance is defined by, 

p T (lcpton) > 20 GeV , y (lepton) < 1 . (7.5) 
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Treatment of V^-boson and &-quark 


olo [fb] 


vnlo (prod.) [fb] 


<?nlo (prod.+dccay) [fb] 


Narrow width, m& = 
Narrow width, m& = 4.7 GeV 
Breit-Wigner, m\, = 4.7 GeV 


12.14(2) 
12.12(2) 
12.08(2) 


19.96(2) 
19.96(2) 
19.88(2) 


20.03(2) 
20.01(2) 
19.95(2) 



Table 4. Cross sections in femtobarns for the s-channel single top process, computed with various levels of 
sophistication. The cuts appropriate for the Higgs search, as described in the text, have been applied. 



In addition wc require that there is at least 20 GeV of missing transverse momentum. Our results 
for the process pp — > t(— > i^e + 6)6 under this set of cuts are shown in Table 4. As anticipated, the 
application of cuts results in a small difference between the two columns of NLO results. We find 
that the NLO cross section including radiation in production and decay stages is approximately 0.3% 
higher than the result for radiation in production only. For similar reasons, the cuts also induce a 
slight reduction in the cross section (by less than 0.5%) when working with the full Breit-Wigncr 
propagator rather than in the narrow width approximation. Wc also see that the effect of including 
the 6-mass is negligible, due to the requirement that at least two jets are reconstructed with transverse 
momenta that are relatively large compared to the 6-quark mass. 

7.3 i-channel single top at the LHC 

The production of single top quarks via the t-channel process is the dominant production mode at 
the LHC. The final state is hard to distinguish from top pair production and a variety of QCD 
backgrounds such as the production of a W boson in association with jets. In order to distinguish 
single top events from the main background processes, a recent CMS analysis [35] relies on accurate 
theoretical predictions for two obscrvables. 

The first is the rapidity of the light jet that is present in the event. Since at LO this light jet 
is produced by the i-channel emission of a W boson from an initial quark, it is found primarily in 
the forward direction, in contrast to the centrally-produced light jets from background processes. The 
second observable, cos 0* is computed by performing a boost to the rest frame of the reconstructed top 
quark and then defining 9* as the angle between the charged lepton and the light jet in that frame. 
For both these obscrvables it is interesting to investigate to what extent the distributions are modified 
by higher order corrections to both the production and decay processes. 

For the study presented here we have adopted the cuts from Ref. [35]. Light and heavy- flavour 
jets are clustered according to the anti-fey algorithm with a distance parameter D = 0.5 and only 
satisfy a fairly loose set of cuts, 

pr(jet) > 30 GeV , y(jet) < 5 . (7.6) 

The cuts on the charged lepton are, 

p T (lcpton) > 20 GeV , y (lepton) < 2.1 , (7.7) 

and no requirement is made on the missing transverse momentum. Following the CMS analysis, we 
demand that exactly two jets are present after jet clustering: the appropriately-charged bottom jet 
and a light (non-tagged) jet. We reconstruct the top quark by combining the leptonically decaying W 
boson (reconstructed perfectly, by assumption) with the 6-jet. If the invariant mass of the W-6-light- 
jet system is closer to mt than the invariant mass of the W-b system then the radiation is assumed to 
occur in the decay and the event is dropped. 
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t channel single top, 7 TeV LHC t channel single top, 7 TeV LHC 




cose* |y(light jet)| 



Figure 8. Distributions of cos 8* (defined in the text) and |yn g ht at the 7 TeV LHC, computed using the set 
of cuts described in the text and adapted from Ref. [35]. 



Our results are shown in Figure 8. We observe that the effect of the NLO corrections in the 
production is a slight change in the shape of both distributions. However, the effect of additionally 
including the radiation in the decay is not significant. 

7.4 Top pair phenomenology 

The phenomenology of top pair production is very rich in comparison with the single top processes 
due to the much higher production cross section. We therefore perform a more detailed study of top 
pair production in this section. 

7.4.1 Top production at the LHC 

We first provide predictions for some basic kinematic distributions that have recently been analyzed 
by the CMS collaboration [36]. The cuts employed in our study are the ones presented in that analysis. 
Both light and heavy-flavour partons are clustered according to the anti-fcy algorithm with a distance 
parameter D = 0.5 and the resulting jets must satisfy the acceptance cuts, 

prQet) > 30 GeV , y(jet) < 2.4 . (7.8) 

We require that two 6-jets are found by the algorithm, out of a total of at least two jets (for the dilepton 
case, with the top and anti-top decay leptonically) and at least four jets (in the lepton+jets case, with 
one top decaying hadronically). The cuts on the charged lepton depend on the decay channel: 

dilepton : p T (lcpton) > 20 GeV , y(lepton) < 2.4 , 
lepton+jets : p T (lepton) > 30 GeV , j/ (lepton) < 2.1 , (7.9) 

but no requirement is made on the missing transverse momentum. 

Wc will also present distributions for the W bosons, top and anti-top quarks that are reconstructed 
according to the following algorithm. For simplicity wc assume that a leptonically decaying W is 
perfectly reconstructed. For a hadronic W decay we consider as candidates all pairs of light jets and, 
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if there are three light jets, also the system of all three together (to account for radiation in the decay). 
The system whose invariant mass is closest to mw is assigned as the hadronic W decay. For the top 
and anti-top quarks, we consider both the system consisting of the W and the appropriate bottom 
quark, and the system that also contains the remaining light jet (if it was not already assigned to the 
W decay) . The combination of assignments that results in invariant masses closest to for both top 
and anti-top quarks is considered the proper solution. We note that, since we produce both top and 
anti-top quarks on-shell, this reconstruction is often perfect. However, in the case that initial state 
radiation in the production stage is merged with a bottom quark in the final state, one of the top quark 
masses may be reconstructed far from its mass shell. We have not tried to remove such configurations 
in the results presented here, although the code is flexible enough to pursue such approaches. 

In Figure 9 we present a selection of observables for the dilepton process. Focusing first on 
the leptonic observables we see that the NLO corrections have a considerable impact on the shape 
of the pt distribution of the charged leptons but that their rapidity distribution is left unchanged. 
Moreover, the effect of including NLO corrections in the decay of the top quark is negligible. Turning 
to the observables obtained after reconstruction of the top and anti-top quark momenta, we see that 
the rapidity distribution of the top quarks is not affected by the QCD corrections. The shape of 
the transverse momentum distribution of a single top quark does receive substantial corrections at 
large pr ~ 250 GeV. In addition, the two NLO curves also differ: the prediction when including QCD 
radiation in the top decay lies between the LO and NLO (production only) curves. 

For the lepton+jets channel we consider the process pp — > t(—> e + iyb)i(—^ e~vb) + X, i.e. we 
consider the leptonic decay of the top quark and the hadronic decay of the anti-top quark. In Figure 10 
we show the invariant mass spectrum of the W + and W~ bosons produced in this process. We note 
that for the W + boson we have made the simplification that the neutrino 4-momentum is known, while 
the W~ boson is reconstructed according to the algorithm above. We see that, at LO, the predicted 
distribution for W~ is identical to W + , since there is no ambiguity in identifying the jets (we assume 
that the &-jets are perfectly tagged). However, at NLO the additional radiation - either before or 
after the top quark decay - may be mistakenly assigned to the W~ boson, leading to the substantial 
change in shape of the distribution. In addition there is a significant difference between the two NLO 
predictions, with radiation in the decay leading to a bigger enhancement in the region rrijj < Mw- 

In Figure 11 we show the invariant mass of the top, anti-top system in both the dilepton and 
lepton+jets channels. This distribution is important for New Physics searches and, in both cases, 
although we observe important NLO effects at high m t j, this is purely due to the treatment of NLO 
effects in the production stage and not in the decay. 

Finally in Figure 12 we present the distributions of two additional observables for the dilepton 
channel - the transverse momentum of the positively charged lepton, pt{^ + ) and the invariant mass of 
the £ + and 6-quark system, m e + b . As first pointed out in Ref. [14], in the context of the LHC operating 
at 10 TeV, the latter distribution is particularly interesting since it exhibits differences between the 
NLO predictions with and without including QCD corrections in the decay (in the region just below 
the LO threshold at m t ). Here we simply note that the same pattern is observed at 7 TeV, under 
typical analysis cuts used for the 2011 run, as was noted in Ref. [14]. In passing we also remark that 
our predictions are in complete agreement with those of Ref. [14] once the appropriate differences in 
input parameters are accounted for. 

7.4.2 The top quark forward-backward asymmetry at the Tevatron 

A feature of top quark production at the Tevatron that has recently received much attention is the 
top quark forward-backward asymmetry [37-40]. The predicted asymmetry is only non-zero at NLO 
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top dilepton channel, 7 TeV LHC top dilepton channel, 7 TeV LHC 




top dilepton channel, 7 TeV LHC top dilepton channel, 7 TeV LHC 




Figure 9. Predictions for various observables in the top pair dilepton channel at 7 TeV LHC, calculated at 
LO (black), NLO in production (blue) and NLO in production and decay (red). The observables are: the pr 
(top left) and rapidity (top right) of the charged leptons and the p t and rapidity of the reconstructed top and 
anti-top quarks. In each case, a single event enters the histograms twice, once each for the particle {I or t) 
and anti-particle (I or F). 

and beyond and, as an example of the utility of our calculation, in this section we provide a prediction 
for the parton-level asymmetry within the fiducial coverage of the CDF detector. 

To that end, we focus on the lepton+jets channel and adapt an analysis approach based on that 
presented in a recent CDF note [40]. Partons are clustered using the anti-fey algorithm with distance 
parameter D = 0.4 and we require that at least four jets are found that satisfy, 

prQet) > 20 GcV , y(jet) < 2 . (7.10) 

In addition, we demand that both the b and b quarks are identified within these jets, with those jets 
detected in the restricted rapidity range, y(b — jet) < 1. For the Ieptonic top decay we require, 



PT(lepton) > 20 GeV , y(lepton) < 1.5 , 



(7.11) 



top lepton+jets channel, 7 TeV LHC top lepton+jets channel, 7 TeV LHC 
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Figure 10. Distribution of the W-boson candidate invariant masses in the lepton+jets channel at the 7 
TeV LHC. The top quark decays leptonically, so the W + candidate (left) is perfectly reconstructed, while 
the anti-top quark decays hadronically leading to ambiguity in the reconstruction of the W~-boson at NLO 
(right). 




Figure 11. Distribution of the invariant mass of the top, anti-top quark system at the 7 TeV LHC, for the 
dilepton channel (left) and the lepton+jets channel (right). 



and at least 20 GeV of missing transverse momentum (from the neutrino). The top and anti-top quarks 
are reconstructed using the same procedure as in the LHC lepton+jets analysis described above. 

As for all the previous results, our predictions are based on the choice of scale [ir = /jlf = m t . 
Ordinarily a NLO prediction is relatively insensitive to the choice of scale but, in this case, since the 
asymmetry is absent in the LO prediction our results depend rather strongly on the choice of scale. 
For this reason we also consider variations of this scale by a factor of two about this central choice, 
i.e. in the range (m t /2, 2m t ), and take the variation in the prediction as an estimate of the theoretical 
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Figure 12. Distribution of the transverse momentum of the positively charged lepton, £ + (left) and the 
invariant mass of the £ + and 6-quark, m e + b (right) in the top pair dilepton channel at the 7 TeV LHC. 



uncertainty due to uncalculated higher orders. For each calculation we compute the asymmetry, 

A = ctnlo (Ay > 0) - °nlo (Ay < 0) 
FB (tnlo (Ay > 0) + (Tnlo (Ay < 0) ' K ' ' 

where the rapidity difference is defined by Ay = y t — yi- Our results are presented in Table 5 where, 
for the sake of comparison, we also include predictions for the asymmetry in the absence of any cuts 
on the decay products of the top quarks. In addition to computing ^4fb inclusively, we also present 
predictions broken down into contributions over various rapidity difference ranges. As can be seen 
from Table 5 the theoretical uncertainties permit excursions from the central value of as much as 60%. 
Had we used the leading order prediction in the denominator of Eq. (7.12) the estimated theoretical 
uncertainties would have been smaller. As such, the use of the NLO result in the denominator is a 
conservative choice. We observe that the asymmetry expected in the fiducial range is smaller than 
that predicted in the full rapidity range, but grows with Ay. The predictions for the two NLO 
calculations, with and without NLO effects in the decay, do not differ greatly. Since the component 
of the calculation that includes radiation in the decay contains no asymmetry one might expect the 
computations that include these effects to result in a smaller value of Ape- In contrast, we find that in 
the region Ay > 0.5 the cross section under the experimental cuts is lower in the presence of radiation 
in the decay, so that the denominator in Eq. (7.12) is reduced and the predicted Afb higher. 

Our results are compared to the CDF results reported in Ref. [40] in Figure 13. 2 The inclusive 
calculation, with no cuts, can be compared in a straightforward manner with the "parton-level" 
corrected results that are presented in the CDF note. In contrast, when the experimental cuts are 
applied, the comparison must be interpreted with considerable caution. Our analysis differs somewhat 
from that used by CDF, notably in the treatment of the top quark reconstruction, but most importantly 
the CDF data is simply background-subtracted and not corrected back to the parton level. In that 
case, our theoretical predictions should be interfaced with a parton shower in order that a full detector 

2 The theoretical predictions in Ref. [40] include a 26% correction to the asymmetry because of elcctroweak contri- 
butions. These corrections will not be included in our results. 
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Rapidity range 


A™ cuts (NLO) 


^™ ts (NLO production) 


,4 F u B ts (NLO prod + decay) 




inclusive 




0fi5+ 028 




045+ 021 

U.U4O_ on 


< 


|Ay| 


< 0.5 


0S6+ 020 
u.uod_ 003 


o m7+ 005 

u - ui ' -0.005 


01 fs+ 009 
u.uio_ 004 


0.5 < 


|Ay| 


< 1 


062+ 013 

U.UUZ_ 028 


0.05218:83 


053+ 025 
u.uoj_ on 


1 < 


|Ay| 


< 1.5 


101+ 060 

U.1U1_ .006 


0.086±g;°?| 


092+ 039 

U.U»Z_ 021 




|Ay| 


> 1.5 


U.iyo_ 036 


o.i42i8;8§S 


14O+0 062 

u.i 2 ±y„ 03 6 



Table 5. Predictions for the top forward-backward asymmetry in the lepton+jets channel at the Tevatron, 
computed without applying any experimental cuts (A P B CUts ) and also when using the cuts described in the 
text (Ap'Jl' 5 ). The uncertainties are obtained by varying the scale in the range (mt/2,2mt). 
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0.5 
0.45 
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Figure 13. The top quark forward-backward asymmetry, Afb, at the Tevatron, after acceptance for correc- 
tions (left) and in the measured fiducial region as detailed in the text (right). Data are taken from CDF [40] 
(Tables IX and XVI) and the MCFM predictions are for inclusive tt production (left) and for the differential 
rate including radiation in both production and decay (right). The central prediction is shown as a solid line 
and the dashed lines represent an estimate of the theoretical uncertainty. 



simulation be performed before making a definitive comparison. Nevertheless, our results can be taken 
as a guide to the level of uncertainty expected in current predictions of the asymmetry as a function 
of Ay. 

We conclude with predictions for an observable that does not depend on reconstructing the top 
quark kinematics and that would not be expected to receive significant corrections due to the addition 
of a parton shower. This observable is the asymmetry observed in the charge- weighted rapidity distri- 
bution of the lepton (i.e. qirji). The lepton inherits the asymmetry of the top quark and its properties 
should be predicted more robustly by our calculation. Using the same procedure for assigning the 
uncertainty as before we find, 

A£g ton (NLO production) = 0.021±g;ggg , A^ ton (NLO prod+decay) = 0.020l°;8o3 ■ ( 7 - 13 ) 
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These predictions should be compared with the most recent CDF result [40] , 

^pg ton (CDF) = 0.065 ± 0.020 . (7.14) 

8 Conclusions 

In this paper we have described in detail the implementation of single top and top pair processes in 
the parton-level integrator MCFM. We have presented an amplitude-level treatment that allows the 
inclusion of the decay at NLO for all processes involving top quarks. We have currently implemented 
this scheme in MCFM for three processes, top-pair production, s-channel single top production, and 
i-channel single top production. 

Our treatment of top-pair production has been rewritten to use this method and to incorporate 
a faster treatment of the one-loop amplitudes. By retaining the mass for the b quark we can work 
in a four-flavour scheme for the i-channel single top process. In this scheme the bottom quark does 
not appear in the initial state, but the production cross section depends logarithmically on the mass 
of the 6-quark. The four-flavour i-channel single top production, implemented using this method and 
presented in this paper is new. 

The features we have described arc complete as of MCFM v6.2. MCFM therefore provides the 
most sophisticated possible NLO treatment of these top production processes within the context of the 
top pole approximation. With these features we can assess the importance of NLO radiative effects in 
production and decay, the effect of the b quark mass, and the effect of an off-shell VF-boson. For the 
distributions that we have examined, these effects turn out to be quite small. 

We look forward to confronting this implementation of NLO theory with data. 
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A Spinor notation 

Our spinor notation is quite standard in the QCD literature, (for a review see refs. [41, 42]). The 
function u±(ki) is a masslcss Wcyl spinor of momentum ki and positive or negative chirality. In terms 
of these solutions of the Dirac equation, the spinor products are defined by, 

(ij) = (i-\j + )=u-(k i )u+(k j ), (A.l) 

[ij] = (i + \n = «+(&)«_(*,) . (a.2) 

We use the convention [ij] = sgn(fc"fc°) (j i}* , so that, 

(ij) [j i] = ■ kj = Sij . (A. 3) 

and 

(ji) = -(ij), \ji] = -[*j]- ( A -4) 

For massless spinors we have the following generalization of the Fierz identity, 

(a\r\b]^ = 2(\a)[b\ + \b](a\) . (A.5) 
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Another useful identity is the Schouten identity, 



(ab) (cd) = (ad) (cb) + (ac) (bd) , 

[a b] [c d] = [a d] [c b] + [a c] [b d] . (A.6) 

We further define, 

(a|y|6] = (a|.p|fr] — > (ap) [pb] for massless p, (A. 7) 

where the decomposition into a pair of spinor products is only valid for a massless momentum p. 

B Calculation of the total width for top decay 

Our calculation retains the NLO corrections to the fully differential top decay rate. This requires 
the inclusion of the correlations with the decay products of the W^-boson as well as the effects of 
extra gluon radiation if present. The calculation of the total width that we detail here sums over the 
polarizations of the W-boson and integrates over the momentum emitted gluon. Although this is not 
a new result, the total width is one of the ingredients in our calculation. In addition, the calculation 
of the one-loop contribution to the differential rate is almost identical to the calculation of the one- 
loop contribution to the total width and our integration over the real radiation counterterm closely 
parallels the integration over the real radiation given below. Our treatment of the corrections to the 
width follows, and in places supplements, the nice discussion given in ref. [20]. 

B.l Phase space for tree graph top decay 

In the rest frame of the top it is straightforward to obtain the two particle phase space for the decay 
t -> W + b, 

d^ 2) {jp t ;pw,Pb) = ( 27r )»-i (2tt)"-i ( 27F )" S " ( pt ~ pw ~ Pb ) 5< $> ~ m ^ ~ m ^ 



(4tt) 2£ 1 
32tt 2 (m 2 t y 



A(l, W 2 ,/3 2 ) d n ~ 2 a w , (B.l) 



where uj = mw/mt, ft = nit,/ nit and A(x,y, z) is defined in Table 1. The matrix element squared for 
this process, summed over spins, is, 

M KA* ( 4. PwPw \ _ 2G F m* 

M,M v y-g +_^j__^=_ / , (B.2) 
with / also defined in Table 1. Using the fact that, 



/ 



taking the limit e — > and including the flux factor (2nit) 1 , we recover the result for the lowest order 
width quoted in Eq. (2.2). 

B.2 Virtual corrections to top decay 

The general structure of the virtual corrections to the process t — > W + b has been outlined in 
section 4.2. Here we will give results for the coefficients Cq, Cq , Cf and Cf 1 that appear in the form 
factor decomposition, Eq. (4.10). 
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Integral function 


Expression 


A (m t ) 

Ao(m b ) 

B (p t ;0,m t ) 

B (p b ;0,m b ) 

B (p w ; m b , b t ) 

C Q {p b ,p t ;0,m b ,m t ) 


ml(\ + \n^ + l) + 0{e) 

m2(I + ln^ + l-ln/3 2 )+0(e) 

7 + 1^ + 2 + 0(6) 

i + ln^ + 2-ln(/3 2 ) + 0( e ) 

± + ln4 + 2 + + 0(e) 

-k^(^ + ln $) +c(/) + H 



Table 6. Table of one-loop integrals obtained from ref. [43]. The functions c"' and b^' are given in the text. 



These coefficients involve the following functions related to the scalar integrals, 



Bo(pi;mi,m 2 ) 



Ca{pi,p 2 ;m ll m 2 ,m 3 ) 



,4-D 



i{2ir) D c r 



id, 



1 



(I 2 — m\ + ie) 



/' 



A-D 



i(2n) D c r 

,4-D 



d°l 



/' 



i(27r) D c r 



d D l 



{I 2 - m{ + ie)((l+pi) 2 - mi. + ie) 

1 



(B.4) 



(I 2 - m\ + ie)((l +pi) 2 - m\ + ie)((l +pi +P2) 2 — w§ + ie) 



where the constant cr has been defined in Eq. (4.11). Explicit expressions for the particular scalar 
integrals that appear in the calculation are given in Table 6. In order to have a compact representation 
for the finite parts of the integrals, we have introduced the following functions, 



uj- 



z u) = Li 2 (l - P_) - Li 2 (l - P+) - Li 2 ( 1 - 



71 



In 2 - In 2 P + 



(B.5) 



The results for the coefficients in the four-dimensional helicity scheme [44] are, 



cf 



- 1 



H + ^ c C/)_4 + ln(/J 2 ) 



LftW [1 - /3 2 ^ 2 - w 2 + /3 4 - 2/3 2 - 6P : 



2P| 

A 
j 3 

p 2 
1 

~p2 

r 3 



P \n(p 2 )~W b (s) 



i(l-^ 2 -/3 2 )^)-/3 2 ln(/3 2 ) 



* ln(/3 2 )[3P 2 + /3 2 w 2 -/? 4 + /3 2 
lf 3 



(B.6) 
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These expressions include the effects of wave function renormalization, as indicated in the upper-left 
diagram of Figure 6. The result in the four-dimensional helicity scheme is, 

,2 

■5 



Z Q = 1 - 9 c r Cf 



-+31n(^ 



(B.7) 



The result for the wave function renormalization is independent of the gauge fixing parameter. 
In terms of these coefficients the matrix element squared is, 



2G F m\ 



nv , PwPW 

y ' 2 



V2 



f + 2g'c r ^ C 



Co/ + (Cf + 0Cf)(£ - 3^ 2 P ) - 6/3a; 2 C fl 



0(4) (B.8) 



Inserting the values for the coefficients from Eq. (B.6) and including the wave-function renormal- 
ization the total virtual result is, 



THVl 



Too 4:C F g 2 c r ( — o 



i [P r p - Ps] + Po [U 2 (l - P_) - Li 2 (l - P+) - Li 2 fl ~ ^ 



y 2 - 2 (y„ + in /?) (y w + r p ) + 2F W in p+ 



2P a 



12w 2 F p P 3 2 



[1 + 4/3 2 - 5/3 4 - w 2 (5 - f3 2 ) + 4w 4 ] P 3 In /3 



(B.9) 



in agreement with Czarnecki [20] after using standard identities between dilogarithms. 
B.3 Real radiation 

The evaluation of the real corrections relies on a factorization of the phase space into a simple form 
that is suited both to the calculation of the total width and the more differential calculations presented 
in this paper. 

B.3.1 Factorization of Phase Space 

For radiation in the decay the appropriate phase space is given by, 
_(3), s d n p w d n p b d n p g 

d * 'to**'**') = ( 27r) »-i ( 27r )n-i (2*)-i 

x 5+(p 2 w - m 2 w )8 + (pl - m 2 b )8 + {p 2 g ){2n) n 5 n {p t -p w -p b -p g ) . (B.10) 

By inserting an integral over a delta function, S n (p t —pw — P)d n P, we can rewrite this in the familiar 
factorized form, 

d^ 3 \p t ;p w ,p b ,p g ) = J — d^ 2 \p t ;p w ,P)d^ 2 \P;p b ,p g ) , (B.ll) 

where P = p b + p g . Working in the center of mass frame of P = p b + p g we have, 
m t (z- ft 2 ) 



Pg 



Pb 



2^ 
ro t Q + /3 2 ) 
2^ 



1, . . . , sin#, cos 6) , 
(1, . . . , — v sin 9, —v cost 



with v 



z-f3 2 
z + f3 2 



mi 



Pt = (E t ,...,0,P t ) , with Et = r7=(l + 2-l/), Pi 



(B.12) 
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Hence P 2 = m 2 z and we can use the result of Eq. (B.l) to write the factorized form, Eq. (B.ll), as, 



d& 3) (p t ;p w ,Pb,Pg) = d 
Zir 



32tt 2 (m 2 t y 



\(1,L0 2 ,Z) 



jn-2 



-i , 



(4tt) 2 



1 



2\l-2e jn-2 



(B.13) 



32tt 2 (m 2 ) e 

By comparing Eq. (B.13) with the lowest order space space in Eq. (B.l) we can write the real emission 
phase space as a leading order normalization multiplying a term containing the real emission degrees 
of freedom, 

d<f> {3) (p t ;pw,Pb,P g ) = d<f> {2) (p t ;p w ,p b ) x [dg(p t ,p w , z)} . (B.14) 
The real emission factor factor [dg(p t ,pw, z)] is given by, 
(m?) 1 "^)* f r X(l,w 2 ,z) 



[dg(pt,pw,z)\ 



64tt 3 



crK) 1 



dz 



A(lV,/3 2 ) 
A(l,w 2 ,z) 



_ «2^-2e rf n-2^ 



z ^(z-/3 2 ) 

-!+£/-„ a 2\l-2e 



dx(l-x 2 Y 



(B.15) 



.A(l, w 2 ,/3 2 ). 

In this equation we have set x = cos0 and cr has been defined in Eq. (4.11). 
B.3.2 Matrix element squared and integration. 

We now present the result for the spin-averaged matrix element squared for the real radiation process, 
t -> W + b + g. The result is, 



au i PwPw \ _ 



f 



2m 2 P Q 



+ f 



V2 

2p b .p g ~ 2p t .Pg 



Pt-P 2 g 



Pb-Pg Pt-Pg Pb-Pg 



(B.16) 



Pt-Pg Pb-Pg 



4UJ 2 (pt.p 2 g+Pb-P 2 g) 2{/3 2 + l)(PfPg- Pb-Pg? 



H Pt-Pg Pb-Pg 



Pt-Pg Pb-Pg 



This expression can be written in terms of the variables of integration in the phase space measure, z 
and x = cos9 (c.f. Eq (B.15)), by using the identities, 



Pb-Pg 



■{Z-(i 2 ) p t .Pg =Eg{E t -P t COS 9) 



(B.17) 



It is convenient to first perform the angular integration over x = cos 9. From the expression for the 
matrix clement in Eq. (B.16) and the form of the phase space in Eq. (B.15), these integrations take 
the form, 



4 ( 

T 



dx (l-x 2 )-*-, 1 — 

l (PfPgY 



(B.18) 



for n = 0, 1 and 2. The results for J , Ji and Ji are given in Table 7. 

Wc shall first focus on the integration of the divergent terms appearing in the first set of parentheses 
[. . .] in Eq. (B.16). After performing the integration over x these terms lead to the following expressions, 
with the integral over z not yet evaluated. For the first term, containing the factor m 2 J (pt.p g ) 2 , we 
are left with, 

2 ^ /VV /" ra , f /z-/3 2 \i- 2 V A(l,w 2 ,z) \!- £ m 2 r n E t „ n . 2 ,i 



g 2 c r C F 



,2 \ £ 



dz 



(z-/3 2 ) 



2\l+2e 



.A(l,w 2 ,/3 2 ) 



1 + 2e-±Y n 
Pi 
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Integral 


Result 


Jo 
Ji 
J2 


1 + 2 

E B Pt{ 


e + 0(e 2 ) 

Y p + e[Li 2 (l-^) 
{l + 2 e fy p + 0( e 2 )' 


+ Y p 2 ]+0(e*)} 

} 



Table 7. Results for the angular integrations specified in Eq. (B.18). 



The expression for the second term, with the factor ml/(pb.p g ) 2 , is, 

,2 \ c fZ m 



g c T Cf 



dz 



g z c T C F I — o 



t / Jf3 

2 \ c r z 





l-2e 4^2 


< A(l,w 2 ,z) \ 




(z - p 2 ) 2 


U(1, W 2 ,/3 2 )J 



ns Jf P (z-^y+^ VA(l,a,2,/3 2 ), 
The corresponding result for the third and final term, proportional to Po/(p t .p g pb-p g ), is, 



S3 = -g l c v C F ( — j 
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2\ 
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A(l,/W) ^ 



z^z-/? 2 )- 1 - 2 ^ 



A(l,o; 2 ,z) 
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1-^1 



0( e 2 )} 
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We may now perform the integration over z for each of these terms. For convenience we give useful 
results for the non-trivial integrands in Table 8. The notation is defined by, 



</(*)>= / dzf(z) 

102 



(B.22) 



Performing the integrals in Eqs. (B.19), (B.20) and (B.21), using the results of Table 8, we obtain 
the following explicit results, 



S x = -4g 2 c r C F 



2 \ e 



2e \pu) °P 3 °P 3 



(B.23) 
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/4P 



S 2 = -4g 2 c T C F ^ + ^ -(1-^)^-^(1-^-^)^ , (B.24) 
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Y P {z) 

z> 

Ps(z) 
z 

Pa(z) 


% (ln(4P 2 ) - %) + |Li 2 (l - P_) - ±Li 2 (l - P + ) 

-Li 2 (l - f^) + ln(P_) ln(l - cu - P_) - ln(P+) ln(P+ - 1 + w) 

±P 3 (1 + /? 2 + 5a; 2 ) - ±u; 2 (2 + W 2 )F W - 
P 3 - u 2 Y w - f3% 

(1+W 2 )y- t7 , P 3 

(1-ou 2 )Y p -uj 2 Y w -P 3 
-u 2 Y w + W P 3 



Table 8. Results for finite z integrals as defined in Eq. (B.22). Notation as in Table 1. 



The remaining finite contributions can be integrated in four dimensions and added to the above results 
to yield the full real contribution to the width, 

a s r? al = T^ACpg 2 ^ ( ^) hf\ - -[P Y p - P 3 ] + P [Li 2 (l - P_) - Li 2 (l - P+) - 3Li 2 (l - 

\mfj Lie L P + 



Y 2 + 2(Y P + ln(/3))(K, + Y p ) + 2 In 



4P| 



Y„ 



W+P 2 

+ 2(1 - p 2 ) [(1 - 1 ) 2 + w 2 (l + P 2 ) - A^]Y W 

+ - [3 - [3 2 + 11/3* - f3 e - 4/3 4 w 2 - 9w 4 + 7^ 2 o; 4 + 6oj 6 ]Y p 

+ 4 [1 - 2/3 2 + /3 4 + w 2 + /3 2 w 2 - 2w 4 ] P 3 InOS) 

+ i [5 - 22/3 2 + 5/3 4 + 9w 2 (l + /3 2 ) - 6w 4 ] P 3 ] ■ 



2 + 21n 



(w)H 



(B.26) 



Adding in the virtual correction as given in Eq. (B.9) we find that the total correction to the width is 
in agreement with Eq. (2.6). 
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